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Single partially confined collapsed polymers are studied in two dimensions. They are described 
by self-avoiding random walks with nearest-neighbour attractions below the 0-point, on the sur- 
face of an infinitely long cylinder. For the simulations we employ the pruned-enriched-Rosenbluth 
method (PERM). The same model had previously been studied for free polymers (infinite lattice, 
no boundaries) and for polymers on finite lattices with periodic boundary conditions. We verify the 
previous estimates of bulk densities, bulk free energies, and surface tensions. We find that the free 
energy of a polymer with fixed length N has, for TV — > oo, a minimum at a finite cylinder radius 
R* which diverges as T — > Tg. Furthermore, the surface tension vanishes roughly as (Tg — T) a for 
T -^T e with a « 1.7. The density in the interior of a globule scales as (Tg - T) with (3 « 0.32. 



Although the behaviour of polymers has been studied 
for many years, there are still a number of open questions. 
At high temperatures or in good solvents repulsive inter- 
actions (the excluded volume effect) dominate the confor- 
mation and the polymer chain tends to swell to a random 
coil. At low temperatures or in poor solvents attractive 
interactions between monomers dominate the conforma- 
tion and the polymer chain tends to collapse and forms a 
compact dense globule. The coil-globule transition point 
is called the O point. 

The behaviour in good solvents is accessible by field 
theory and thus well understood Q, e.g. the partition 
function scales as Z ~ p N N 1 ^ 1 at T > Tg. In poor 
solvents, much less is known exactly. For a long time it 
was e.g. widely believed that the partition sum scales 
in the same way, although already Lifshitz mean-field 
theory J|, |j| had predicted for T < Tg that there should 
be also a surface term, 

Z(N,T) ~ a N b- N "N^- 1 , (1) 

with s = (d—X)/d and b > 1. First numerical indications 
for this were given by ||, |5|, [| , and theoretical arguments 
in favour of Eq. (|l|) were given in Q|. More recently, 
Eq. was checked numerically in M . But it seems that 
there exists still no proof for Eq. (EJ), and it is not known 
how the surface tension (related to b) scales for T — > Tg. 
Moreover, the exponent 7 is not known ^, |^], and it 
is not even clear whether it is universal in the collapsed 
phase. 

In (|] we had studied two different geometries - glob- 
ules on infinite lattices and polymers filling finite lattices 
with periodic boundary conditions - where surface ten- 
sion should make very different contributions. In the 
present work we study a third geometry, namely poly- 
mers on the surface of an infinitely long cylinder, where 
the surface tension should make yet another contribu- 
tion. Our main motivation was originally to have a sim- 
ple geometry with straight "surfaces" so that the surface 
tension and its T-dependence can be more easily mea- 
sured. But we found also another interesting effect: the 
bulk fugacity p depends non-monotonically on the cylin- 
der radius. It has a minimum at a finite radius R* which 
diverges for T — > Tg. Thus a collapsed polymer, when 
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FIG. 1: Schematic drawing of a polymer chain growing on the 
surface of a cylinder. 



wrapped around an elastic cylinder, will tend to squeeze 
it when R> R* and widen it when R < R*. 

A collapsed polymer on the surface of an infinitely 
long cylinder is modeled as an interacting self-avoiding 
random walk (ISAW) on a square lattice with periodic 
boundary conditions in the vertical direction, as shown in 
Fig. [|. The vertical lattice size is h = 2-kR. The collapse 
is mediated by nearest neighbour attraction between 
monomers, described by a Boltzmann factor q — e~@ e 
for each pair of non-bonded monomers occupying near- 
est neighbour sites on the lattice. Here f3 = 1/kT and 
e < is the attractive potential. The partition sum is 
therefore 

Z = E ^ > ( 2 ) 

walks 

where m is the number of non-bonded occupied nearest- 
neighbour pairs. We are interested in polymers which 
form collapsed globules, therefore q > qg with qg = 

As argued in [p], |3) and ver ified in @, §, below the 
point, the polymer forms a globule with uniform 
monomer density p for large N and p — > as T — > Tg. 
Thus, the free energy for a collapsed polymer with chain 
length N in infinite volume consists (except for a loga- 
rithmic term resulting from the last factor in Eq. ([!])) of 
two parts: An extensive bulk contribution cx N, and a 
surface contribution oc N^ d ^ 1 ^/ d . For d = 2, this reads 
(up to a factor 1/T which will be suppressed in the fol- 
lowing) 

-InZjvO) ~ Poo(q)N + a{q)VN + 0{\nN) (3) 

where p^ is the chemical potential per monomer in an in- 
finite chain, a is related to the surface tension (=free en- 
ergy per unit of perimeter length) a and to the monomer 
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density p by 

The latter follows by assuming the globule to be a disk 
with radius R g fixed by N — nR 2 p. 

Surface terms are absent from the free energy if one 
has a finite L x L lattice with periodic boundary condi- 
tions, and if one chooses the polymer length such that the 
entire lattice is filled with finite density. More precisely, 
consider the quantity ln^Ar(g) + pN . For T < Tg this is 
non-convex, and for p = pl (q) it has two peaks of equal 
height, corresponding to the coexistence of two phases: 
One phase with very short chains and density 0, and 
another phase with density plUi)- For L — > oo, pl(q) 
and Pl{q) tend monotonically || to Poo(q) and p. 

Finally, on the surface of an infinitely long cylinder 
with perimeter h, a collapsed polymer forms a cylindrical 
blob with length L = p~ 1 N/h and surface free energy 
2ah, 

-\nZ N (q,h)^p h {q)N + 2a{q)h. (5) 

Notice that the chemical potential per monomer will in 
general depend on h, even in the limit N — > oo. For 
h — > oo we expect again to obtain the bulk limit, 

Hh{q) ->■ Moo(e) h->oo. (6) 

In principle we should also expect cr{q) and p(q) to de- 
pend on h. But we shall not study the detailed behaviour 
of p, and we shall neglect the /i-dependence of a(q), since 
the surface term is a small correction anyhow. We also 
might expect logarithmic terms in N to be present (just 
as also in the previous geometries), but they will also be 
neglected. 

Our aim is to verify Eq. (|5|) numerically, to estimate 
the functions pn{q) and and to see whether they are 
compatible with the results for the other two geometries. 

The pruned-enriched-Rosenbluth method (PERM) has 
been employed to study collapsed 2-d polymers in |^, [l0| . 
In Q the predicted behaviour was verified numerically 
for globules on infinite lattices and for polymers filling 
finite lattices with periodic boundary conditions. More 
precisely, finite lattices were studied first, and from them 
estimates were obtained for /Xoo and p by extrapolation. 
Given these, subsequent simulations on infinite lattices 
gave rather precise estimates for a(q). But the latter was 
not done systematically for a wide range of temperatures, 
in particular the behaviour of cr{q) for q — > qg was not 
studied. Thus we want to obtain also the latter in the 
present paper from cylinder geometry simulations. 

For the first simulations, we choose q = 2.4 which is 
deep inside the collapsed regime. From Eq. (||) we expect 
that plots of In Zx{q, h) + pN against N should give hor- 
izontal lines for N — > oo, if and only if p = pu{q). This 
gives estimates of Ph{q)- Using them, we plot in Fig. || 
hxZw(q,h) + ph(q)N versus N for various values of the 
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FIG. 2: \nZ N {q,h) + p h (q)N versus N for q = 2.4 and for 
various values of the cylinder circumference h. The values of 
p were determined by demanding that these curves become 
horizontal as N — > oo. The peaks seen in the curves for h — 16 
and 20 are statistical fluctuations. 
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FIG. 3: ph(q = 2.4) as obtained from Fig. ^| plotted against 
h. There is a minimum at h — 6.13±0.01 which is determined 
by fitting a parabola (dashed curve) going through the three 
lowest points (which are at h = 5, 6, and 7). The dotted line 
is the asymptote poo(q) = —1.3211(1) approached for h —* oo. 



cylinder circumference h. For small h we find perfect 
agreement with our expectations. For larger h we see 
that PERM has problems in sampling long chains cor- 
rectly, as should be expected. Fluctuations similar to 
those seen in the curve for h = 16 prevent us from go- 
ing to much larger h, except for q close to qg where the 
efficiency of PERM increases. 

Values of Ph{q), again for q = 2.4, are plotted in Fig. ||. 
For large h, we obtain Poo(q) — —1.3211(1) in perfect 
agreement with |^]. But, surprisingly, the dependence of 
Ph{q) on h is not monotonic. Instead, we see a minimum 
at a finite value of h. In our simulations h is an integer, 
but we can interpolate to continuous values by fitting a 
parabola, ph(q) — Pmm + b(h — h*) 2 , through the three 
lowest points, as shown in Fig. 0. Such a fit gives h*(q — 
2.4) = 6.13(1) and ^ min (q = 2.4) = -1.3243(1). 
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FIG. 4: A{q,h) versus 2h for g = 2.0 and 2.4. The dashed 
line has slope -a — —0.105. The points for h = 16 and h — 20 
have large errors for g = 2.4, as also seen from Fig. ^. Results 
for q = 2.0 (insert) show the difficulty of measuring surface 
tension near 0-point. 



Finally, in order to obtain the surface tension, we plot 
in Fig. [| 



A(q, h) = lim [lnZ/v(g, h) 

N—>oo 



(7) 



against 2h. The limit is of course obtained by using 
the parts of the curves in Fig. ^ where they are hori- 
zontal. From Eq. (@) we expect this to give a straight 
line with slope —a. The result shows that a = 0.105(5) 
at q = 2.40. Actually, the line is slightly curved, show- 
ing that systematic errors are not negligible. For large h 
they seem to arise from insufficient sampling. For small 
h they result from the neglected terms discussed above. 
The present estimate should be compared to the esti- 
mate a = 0.119(5) from ||. Although the discrepancy 
is about 2 standard deviations, we consider this as rea- 
sonable agreement. Results for q = 2.0 are also shown 
in Fig. [I], and illustrate the large corrections near the 
0-point. 

Analogous simulations as shown in Fig. were also 
done at other values of q. The lengths N depended 
strongly on q. Close to the 0-point the algorithm is very 
efficient and we obtained reliable high statistics data for 
N up to 15,000. The more the polymer is collapsed, 
the more difficult are the simulations. We stopped at 
q = 2.4. As seen from Fig. 2, the errors decrease when 
h becomes small, most likely because the walk is effec- 
tively Markovian with memory ~ h 2 , and PERM is most 
efficient when memory is short. In all those simulations 
the results were similar. fJ.h(q) bad a minimum at a fi- 
nite h*(q), and A(q, h) versus h was roughly linear with 
a negative slope. Results of h*{q), /x m in(g)) p(q) and cr(q) 
for a wide range of q are shown in Fig. ^-|| 

In Fig. H we can see that h* (q) increases as q decreases. 
This is to be expected since the transition is contin- 
uous, with a diverging correlation length for T —* Tg. 
Apart from the lattice constant (which is set to 1 here), 
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At h — h*(q), the chemical potential is minimal, 
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FIG. 6: jUmin and /ioc plotted against (q—qe) for various values 
of q. The slope of the straight line is roughly equal —0.34. /Xoo 
were obtained from extrapolation to large h. 



the correlation length is the only relevant length scale. 
The divergence of h*(q) when approaching the point 
thus indicates that the existence of the minimum is not 
a lattice artifact and is presumably a genuine feature of 
the continuum theory. The simulations are not precise 
enough to allow a meaningful fit to a power law and a de- 
termination of a critical exponent. This exponent would 
also be effected strongly by the uncertainty of Tg. To 
reduce the latter to a minimum, we performed extensive 
simulations of free (i.e. infinite lattice) ISAWs near the 
estimated point, with the result qg = 1.9487(5) which 
is in good agreement with the best previous estimate [jl0| . 

In Fig. | we plot ^ min (g) and ^oo(g) versus (q-qg). We 
see that both decrease linearly with temperature. The 
difference between them depends weakly on the temper- 
ature. It seems to vanish for T — > Tg. 

The monomer density p inside the "globule" is con- 
stant for N 3> h 2 3> 1 and fixed temperature. It can be 
measured by measuring the ^-coordinate (parallel to the 
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FIG. 7: Log-log plot of p, as obtained from the growth of 
(|x]v|) with N for large N and h, against (g — qg). For this 
plot we used qg = 1.9487. The dashed line is p oc (g — ge) 0,32 . 



cylinder) of the end point of the chain. Since the chain 
has to grow either into the positive or into the negative 
direction, we have 



d(\x N \) 
dN 



(hp)' 



N 



(8) 



For small values of h there are finite size corrections to 
the density, but it is easy to extrapolate to h — > oo. The 
resulting densities are plotted versus q in Fig. They 
show a power law 



p ~ (q- q g f 



(9) 



with f3 = 0.32 ± 0.02. We are not aware of previous 
estimates of (3, either theoretically or by simulations. 

The values of the surface tension a(q) obtained by the 
same way as shown in Fig. ^ are shown in Fig. ^. They 
follow a power law, a ~ (q — qg) a , with a rj 1.7 ± 0.1. 
For the parameter a this implies that 



(« - qe) a - p/2 ~ (g - qe) 



1.54±0.1 



(10) 



This is incompatible with the theoretical prediction a 
(Tg — T) 7 / 6 quoted in §. 



In summary, we have studied collapsed polymers on 
the surface of a cylinder. Using the PERM algorithm we 
could simulate chains of lengths between 3000 (far below 
Tg) and 10 4 (near Tg) for a wide range of temperatures. 
We have verified that the concept of a surface tension 
applies to 2D collapsed polymers, and have determined 
its scaling near the point. We also determined the 
scaling of the monomer density in the collapsed phase, 
as T — > Tg. In addition, we found a surprising non- 
monotonic dependence of the chemical potential on the 
perimeter of the cylinder. It shows that for any temper- 
ature there is a special value of the perimeter where the 
free energy of the polymer is minimal. The fact that this 
value diverges when the O-temperature is approached 
suggests that this is not an artifact of the square lattice 
but a generic feature. 



0.1 



0.01 



0.001 



0.1 



(q-qe) 



FIG. 8: Log-log plot of a, as obtained from Fig. ^j, against 
(q — qe)- For this plot we used qg — 1.9487. The dashed line 
is a oc (q - qg) 1 ' 7 . 
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